Voltage Control of Exchange Coupling in Phosphorus Doped Silicon 
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Motivated by applications to quantum computer architectures we study the change in the exchange 
interaction between neighbouring phosphorus donor electrons in silicon due to the application of 
voltage biases to surface control electrodes. These voltage biases create electro-static fields within 
the crystal substrate, perturbing the states of the donor electrons and thus altering the strength 
of the exchange interaction between them. We find that control gates of this kind can be used to 
either enhance, or reduce the strength of the interaction, by an amount that depends both on the 
magnitude and orientation of the donor separation. 

PACS numbers: 03.67.Lx, 71.55.Cn, 85.30.De 
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I. INTRODUCTION 

Phosphorus donors in silicon have been the subject of 
increased research interest in recent years due to their 
status as either nuclear P| , or electron [^i] -spin qubits in 
various proposals for a scalable quantum computer. In 
both the cases the exchange energy between neighbour- 
ing donor electrons is of fundamental importance as the 
mediator of the the qubit coupling. Although proposals 
exist for implementing quantum information processing 
with this exchange interaction fixed Q, it is none the 
less desirable to have some control over this parameter, 
particularly for the readout process. Recent studies have 
shown that the magnitude of the exchange interaction 
for donor electrons in silicon is strongly dependent on 
not only the magnitude of the donor separation, but also 
on the relative orientation of the donors within the silicon 
crystal lattice 

In a previous article |^ the effect of J-gate biases on the 
exchange coupling was computed for donor separations 
along the [100] axis only. Here we extend the calcula- 
tions to include separations in other orientations relative 
to the host silicon lattice, in particular separations along 
the [110] and [111] crystallographic axes. A two-qubit 
Kane device, Fig^ consists of two phosphorus donors 
at substitutional sites in a silicon substrate at a depth 
20nm below a 50nm layer of silicon oxide. On top of 
this oxide, and between the donors is a metallic J-gate 
electrode, which, for the purposes of these calculations is 
assumed to be 7nm wide, and infinitely long. A grounded 
back plane lies at a depth of 600Abelow the silicon-oxide 
interface. The distance R, between the donors is var- 
ied as is the orientations relative to the silicon crystallo- 
graphic axis. The electro- static potential created inside 
the device by the application of a voltage bias to the J- 
gate electrode, while the remaining electrodes are held at 
ground, is calculated by means of a commercial package 
that solves the Poisson equation for such semiconductor 



systems 0, details of these calculations for Kane type 
devices can be found in ref m. 
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FIG. 1: The Kane architecture based on buried phosphorus 
dopants in a silicon substrate 



II. DONOR ELECTRON WAVE-FUNCTIONS 

In the absence of a perturbing potential the electron 
wave-function for a phosphorus donor centred at a posi- 
tion R/2, in a silicon substrate can be expressed in the 
Kohn-Luttinger [9|,^3 form: 

iI>{y - R/2) = Fl'°'%r ~ R/2)e'^^-('-^)M^(r). (1) 



The functions ■u^(r) are the periodic part of the Bloch 
functions for the pure silicon crystal, with wave-vectors 
located at each of the degenerate conduction band 
minima. The envelope functions , 



(2) 
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are correctly normalised non-isotropic hydrogen like 
wave- functions, with effective Bohr radius a±. The above 
example is for functions localised around the conduction 
band minima k±z = 27r/a(0, 0, ±0.85), with a = 5.43A 
the lattice spacing for silicon. The parameter 7 = aj^/a||, 
where the values a±_ = 25.09A, a|| = 14.43A, are trans- 
verse and longitudinal effective Bohr radii respectively, 
and are determined variationally 0, . The subscripts 
refer to the electronic, orbital and magnetic {n,l,m) 
quantum numbers respectively. 

To control the strength of the exchange coupling it is 
proposed that a voltage bias be applied to a surface "J- 
gate" electrode, see FigQ] This voltage bias produces 
an electric potential within the substrate, perturbing the 
donor electron wave- functions, and thus altering the ex- 
change coupling. To model this process we have calcu- 
lated the electro-static potential produced inside the de- 
vice by the application of a J-gate potential. The poten- 
tial matrix is then calculated in the basis of generalised 
Kohn-Luttinger states given in Eq^ using hydrogenic 
envelope functions up to and including n = 7, a basis of 
140 states in total. The Hamiltonian is then diagonalised 
in this basis to find the ground-state donor electron wave- 
function: 

1/^(r - R/2; V^) ^ ^ Cn,Lrr.{V) (3) 

71.1,771 

X J2 - R/2)e''^-('-^'/^u^(r). 

In the Heitler-London approximation the two donor 
system is treated as a correctly symmetrised product of 
single donor wave- functions. The fermion anti-symmetry 
of the overall wave-function ensures that the spin sin- 
glet (s), and triplet (t) states are respectively represented 
as even and odd superpositions of single electron wave- 
functions. This approximation is valid when the donor 
separation is large compared to the effective Bohr radii 
of the single electron wave- functions. For donors centred 
at positions ±R, the spatial part of the triplet/singlet 
wave- functions are thus: 

^t/,ir,,V2;V) = -i={^(ri-R/2;T/)V'(r2+R/2;F) 

± V(ri+R/2;F)V'(r2-R/2;T/)}. 

(4) 

III. EXCHANGE COUPLING 

In the Kane proposal for a phosphorus nuclear spin 
quantum computer in silicon 1], the effective Hamilto- 
nian between neighbouring donors is given by 

Heff = g„/x„B«i -I- cr.f,2) + 3AisS(crfi + a'-^) (5) 

Here the a represent the usual Pauli operators, and the 
subscripts n, e represent operation on nuclei and electrons 



respectively. The controllable parameters, through which 
quantum information processing is implemented, are the 
contact hyperfine coupling A, between the donor nucleus 
and its associated electron, and the exchange coupling J 
between neighbouring donor electrons. Both these quan- 
tities are to be controlled through manipulation of the 
electron wavefunction via the application of bias volt- 
ages to control gates. The contact hyperfine coupling is 
dependent on jV^CO)!^, that is, the probability density of 
the electron wave-function located at the position of the 
nucleus. This can be changed by drawing the electron 
wave-function away from the nucleus via the application 
of a positive bias to a A-gate pT^, T^. The exchange 
coupling, as can be seen from Eq^l is equal to a quarter 
of the energy difference between the two electron spin 
single, and spin triplet states. In the Heitler-London ap- 
proach, this is evaluated by calculating the energy differ- 
ence between the even and odd superpositions Eq^ The 
origin of this difference in energy between the spin singlet 
and triplet sates is due to the fact that the even superpo- 
sition state has a higher probability density in the region 
located between the two donors, a region that has a rel- 
atively low potential. The odd superposition state has 
a lower probability density in this region which is com- 
pensated by a higher probability density in the higher 
potential regions located outside the two donors. This 
energy difference can be manipulated by application of a 
voltage bias to a J-gate located above and between the 
two donors. Application of a positive bias will decrease 
the potential in the region between the two donors, thus 
increasing the exchange coupling while a negative bias 
will have the opposite effect. 



IV. RESULTS 

It has been well documented that in the absence of a 
bias potential, the exchange coupling between neighbour- 
ing phosphorus donor electrons in a silicon substrate is 
strongly dependent on both the magnitude and the ori- 
entation of the donor separation with respect to the host 
silicon lattice 0, |^ . 

With this in mind it is not unexpected that this orien- 
tation dependence remain in the presence of a bias po- 
tential. Indeed we find that the presence of this potential 
enhances the oscillations observed in the zero bias case. 
This is illustrated in Fig [3 for separations of increasing 
magnitude along three high symmetry crystallographic 
axes. 

It is evident from these plots, that application of a 
voltage bias to a J-gate can exert significant control over 
the strength of the exchange coupling, both increasing 
and decreasing the strength depending on the polarity of 
the bias. The range over which the exchange coupling can 
be changed is however extremely sensitive to the donor 
separation. This is well illustrated in Fig|31where we have 
plotted the value of the exchange coupling as a function 
of applied voltage bias for donor separations of different 
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FIG. 2: Exchange coupling as a function of donor separation 
along high symmetry crystal directions for different J-gate 
biases. The points represent the position of fee substitutional 
sites. 



magnitudes along the same crystallographic axes. 

For all separations studied it was found to be possible, 
via the application of a suitable negative voltage bias, 
to reduce the exchange coupling to zero, and thus to 
turn the inter-qubit interaction off. Indeed calculations 
at large negative bias, particularly for large donor sepa- 
rations, give negative values of the exchange coupling, a 
result that is unphysical for two electron systems. We at- 
tribute this result to a breakdown of the Heitler-London 
approximation in this regime. 

The speed of two-qubit gates, for both electron and nu- 
clear spin architectures increases with the strength of the 
exchange coupling and so it would be desirable to achieve 
the largest exchange coupling possible. To this end it is 
clear that a donor separation along the [100] axis would 
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FIG. 3: Semi-Log plot of the exchange coupling as a func- 
tion of applied gate bias along high symmetry crystal axes 
for different separation magnitudes. Because of the discrete 
positioning of possible substitutional sites, the exact mag- 
nitude of the separation is different for each direction with 
R = 10.86nm in the [100] direction, R = ll.Slnm for the 
[110] separation and R — 11.28 for the donor separations 
along the [111] axis. 



be preferable, with the magnitude of the separation being 
as small as possible. Of course there are limits on the pre- 
cision to which these donors can be placed, and on the 
separations that can be achieved. Technical issues will 
probably set a lower bound on reliable donor separation 
of approximately 2Gnm. Due to the process by which the 
silicon is doped in the so called "Bottom up" approach 
|13^. it is possible to place donors with atomic precision 
in the same [001] plane (that is a plane perpendicular to 
a crystallographic [001] axis) reliably using atomic force 
microscopy techniques and thermal incorporation. Dur- 
ing the process of silicon over growth however the phos- 
phorus atoms are likely to diffuse within this plane by 
up to several nanometers, diffusion out of the plane is 
expected to be significantly less. 

With this in mind we have calculated the expected 
exchange coupling for donors located at various substi- 
tutional sites within the [001] plane, with a separation 
magnitude of approximately 20nm for different donor ori- 
entations, the results are presented in Fig|Sl The first 
plot shows the exchange coupling for donors separated 
by a vector Rg = i?(cos 9, sin 9, 0), that is the donors are 
in the same [001] plane, and when 9 = are separated 
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R ~ 20nm 




FIG. 4: Schematic of the crystallographic structure of a [001] 
plane in sihcon, showing donors separated by vectors along a 
[100] axis (Rioo), and a [110] axis (Rno). For two donors in 
the same [001] plane, we can define the separation Re relative 
to the [100] axis. 



by vector in the [100] direction. The second and third 
plot shows the strength of the exchange couphng when 
the two donors are not in the same [001] plane. In the 
second plot one of the donors has been displaced to its 
nearest neighbour substitutional site and the donor sepa- 
ration is now R = Rg + a/4(l, 1, 1). In the third plot the 
donor is displaced by one lattice constant perpendicular 
to the [001] plane, R = Re + a(0, 0, 1). The plots show 
that even a small displacements from the ideal [100] sepa- 
ration significantly decrease the strength of the exchange 
coupling, indeed a displacement of just one lattice site 
leads to a reduction of approximately one half. 



CONCLUSIONS 



We find that strength of the exchange coupling be- 
tween neighbouring phosphorus donor electrons in a sili- 
con substrate can be significantly increased or decreased 
by the application of a bias voltage to a surface "J-gate" 
electrode placed above and between the donor qubits. 
The bare (zero-bias) coupling is strongly dependent on 
the magnitude and orientation of the donor separation, 
and this dependence is amplified by the application of the 
voltage bias which has the affect of altering the coupling 
strength by a voltage dependent factor. Regardless of the 
donor orientation the exchange coupling strength can be 
reduced to zero, via the application of a negative J-gate 
bias, and so ideal donor placement is determined by the 
maximum value to which the coupling can be increased. 
To maximise the achievable coupling strength the donors 
would ideally be separated along a crystallographic [100] 
axis, however even small deviations from this ideal place- 
ment, of the order of one lattice spacing, can decrease the 
coupling strength by more than a factor of two. 
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FIG. 5: Exchange coupling as a function of donor separation 
angle from the [100] direction for donors separated by ap- 
proximately 200 A in various planes. In the first figure both 
donors are in the same [100] plane, and the angle 9 refers 
to the donor orientation with respect to separation along the 
[100] axis, as shown in Fig|l] In the second figure one donor 
has been moved to a nearest-neighbour site, a displacement 
of (5 = a/4(l, 1, 1) relative to the positions of the first plot. 
The final plot contains data for which the second donor has 
been displaced by one lattice constant in a direction out of the 
plane, 5 = a(0, 0, 1). Here the points represent the position 
of substitutional fee sites. The exact magnitude of the donor 
separation at these points varies slightly due to the discrete 
distribution of such sites. 
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